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ABSTRACT 

Characterization  of  3-D  underwater  light  fields  from  above  the  sea  surface  requires  passive  and  active  remote  sensing 
measurements.  In  this  work,  we  suggest  the  use  of  passive  ocean  color  sensors  and  lidar  (Light  Detection  and  Ranging) 
to  examine  the  vertical  structure  of  optical  properties  in  marine  waters  of  the  Northern  Part  of  the  Gulf  of  Alaska 
(NGOA).  We  collected  simultaneous  airborne  remote  sensing  reflectance  (Rrs)  in  the  spectral  range  443-780  nm 
(MicroSAS,  Satlantic)  and  lidar-derived  volume  backscattering  (J3)  profiles  (0-20  m  depth,  wavelength  =  532  nm)  during 
August  17  2002  in  shelf  waters  situated  south  of  Kodiak  Island  off  Alaska  (57.48°-58.04°  N,  152.91°- 151.67°  W).  We 
evaluated  the  spectral  response  of  Rrs  to  perturbations  on  vertical  distribution  of  p  by  comparing  the  spatial  variability 
between  aggregated  (250  m  horizontal  resolution  x  1  m  vertical  resolution)  spectral  ratios  and  different  lidar  statistics 
per  bin  (Maximum  ft  per  bin,  mean  ft  per  bin,  Pm,  standard  deviation  of  ft  per  bin,  /^td>  integrated  p  per  bin,  $nt)  or 
group  of  bins  (lidar  volume  extinction  coefficient  of  p  between  0  and  5  m  depth).  Sub-surface  changes  of  pm,  pm t,  and 
Ps td  were  mainly  correlated  with  Rrs  (490)/Rrs  (555)  variability  along  the  flight-track  (Semi-partial  correlation 
coefficients  =  0.12  to  0.21).  Our  results  evidenced  linkages  between  above  and  below-sea  surface  optical  properties  that 
can  be  used  to  derive  water  optical  constituents  as  a  function  of  depth  based  on  combined  passive-active  data. 

Keywords:  lidar,  ocean  color  sensors,  water  visibility,  visible  spectrum,  active  sensors,  remote  sensing,  vertical 
structure,  underwater  light  field  models 

1.  INTRODUCTION 

Passive  and  active  optical  remote  sensing  systems  have  inherent  limitations  for  reconstructing  optical  light  fields  within 
the  upper  oceanic  layers  (-0-30  m  depth).  Inherent  and  apparent  optical  properties  in  the  first  optical  depth  (-20  m  = 
1  !Kd,  where  Kd  is  the  vertically  diffuse  attenuation  coefficient  for  downwelling  light)  can  be  derived  from  inversion  of 
above-water  remote-sensing  reflectance  (Rrs  (X,  0+))  (spectral  range  =  400-700  nm)  measurements  obtained  by  passive 
spectrometers1  If  time  and  geographic  location  are  unchanged,  variability  of  Rrs  (X,  0+)  is  related  to  concentration  of 
different  optically-active  compounds  (particulate  and  dissolved)  and  vertical  distribution  of  those  components7,  .  In  other 
words,  ocean  color  observations  from  airborne  or  satellite  passive  sensors  offer  a  vertically-integrated  view  of  optical 
constituents  trough  the  water  column  and  may  not  be  able  to  discriminate  depth-related  optical  features  such  as  thin 
layers4. 

In  contrast  to  passive  remote  sensing  technologies,  active  sensors  such  as  lidar  (Light  detection  and  Ranging),  can  obtain 
optical  measurements  deeper  than  one  optical  depth  and  can  resolve  vertical  differences  on  signal  strength  by  using  short 
laser  pulses  coupled  with  high  speed  time-gated  detectors'.  Typical  lidar  systems  used  from  airplanes  (e.g.,  FLOE, 
fluorescence)6  and  satellites  (e.g.  CALIPSO)8  have  fewer  wavelengths  compared  to  ocean  color  passive  sensors,  one  or 
two  in  the  visible  (e.g.,  excitation  wavelength  in  FLUOR  is  432  nm)  and  one  in  the  IR  (e.g.,  1 100  nm)  spectrum.  This 
spectral  paucity  represents  a  constraint  to  discriminate  targets  with  distinct  absorption  and  scattering  signatures.  Also, 
differentiation  of  optical  components  is  challenging  when  lidar  measurements  (e.g.,  FLOE)  are  based  on  backscattering 
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since  the  relatively  wide  dynamic  range  of  backscattenng  strength,  and  consequently  possible  optical  identities6.  The 
inclusion  of  polarizers  prior  to  detection  has  been  an  alternative  to  minimize  the  spectral  limitations  of  some  lidar  sensors 
(e.g.,  analysis  of  depolarization  ratios  using  CALIPSO  profiles)9. 

In  the  present  study  we  evaluate  the  potential  use  of  concurrent  passive  and  active  optical  measurements  to  retrieve  the 
vertical  profile  ‘shape’  of  inherent  optical  properties  (lOPs)  in  shelf  waters  of  the  northern  part  of  the  Gulf  of  Alaska 
(NGOA).  Spatial  patterns  of  ocean  color  and  lidar  measurements  are  expected  to  be  linked  due  to  changes  on 
backscattcred  photons  at  green  wavelengths.  We  present  preliminary  results  showing  how  different  lidar-backscattering 
parameters  are  related  to  Rrs(A,,  0+)  spectral  variability  based  on  aerial  surveys  performed  during  summer  and  when 
oceanographic  conditions  favor  the  formation  of  diverse  planktonic  layers  in  the  euphotic  zone  (0-50  m  depth). 

2.  METHODS 

2.1  Aerial  surveys  and  flight  mission  settings 

Airborne  spectral  upwelling  radiances  ( Lu )  and  downwclling  irradiance  (Ed)  in  the  visible  spectral  range  (400-700  nm), 
and  lidar  backscattering  (/?)  data  (green  laser  at  532  nm)  were  gathered  over  waters  of  the  eastern  shelf  of 
Afgonak/Kodiak  Islands  (57.48°-58.04°  N,  152.91°-151.67°  W)  during  August  17  of  2002  (Fig.  1). 


1 53°W  481  361  24'  12'  152°W 

Longitude 


Fig.  1.  Aerial  survey  over  the  NGOA  shelf.  Start  (s)  and  end  (e)  locations  during  August  17,  2002  flight  mission  are 
indicated.  KI:  Kodiak  Island,  NGOA:  Northern  Part  of  the  Gulf  of  the  Alaska.  The  full  dataset  of  remote  sensing 
measurements  was  obtained  between  12:54  and  13:22  pm  local  time  and  over  waters  deeper  (>50  m  depth)  than 
penetration  depth  of  the  lidar  system. 

Optimal  flight  weather  conditions  (i.c.,  cloud-free  skies,  wind  speed  <  4  m  s'1)  were  checked  a  priori  to  maximize  the 
number  of  comparisons  between  passive  and  active  optical  measurements.  For  the  whole  aerial  survey,  the  flying  altitude 
and  speed  was  standardized  at  305  m  and  247  km  h'1,  respectively.  Based  on  these  average  flight  characteristics,  we 
collected  104  passive  radiometric  measurements  (upwelling  and  downwelling)  and  5.6  104  lidar  shots  along  a  total 
distance  of  108  km  and  during  28  minutes.  MicroSAS  and  lidar  data  were  geo-locatcd  every  1  minute  during  the  full 
survey. 
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2.2  Passive  optical  measurements 

2.2.1  Optical  sensors 

Measurements  of  Lu  and  Ed  were  performed  at  41 1,  443,  491,  509,  553,  665,  and  780  nm  (10  nm  bandwidth)  with  a 
spectrometric  Micro  Surface  Acquisition  System  (MicroSAS).  These  wavelengths  were  designed  to  match  spacebome 
ocean  color  sensors  such  as  SeaWiFS.  MicroSAS  has  two  digital  optical  sensors  (L:  OCR-507-R03A,  Ed\  OCR-507- 
ICSA,  Satlantic  inc.,  Canada).  As  specified  by  the  manufacturer,  field-of-view  of  the  upwelling  radiance  sensor  is  28°  in 
air,  and  has  a  typical  saturation  of  5  pWcm^  nm4.  Based  on  this  FOV  and  a  sampling  rate  of  360  observations  per 
minute,  typical  pixel  size  of  Lu  MicroSAS  measurements  was  11.5  m  (along-track)  by  200  m  (across-track). 
Downwelling  irradiance  data  were  used  to  discriminate  cloudy  patches  and  calculate  Rrs(^,  0+)  in  each  location.  The 
irradiance  sensor  has  a  typical  saturation  of  300  pWcm  2  nm4  and  a  noise  equivalent  of  2.5 1 0^  pWcm  2  nm 

2.2.2  Atmospheric  corrections 

A  quasi-singlc-scattcring  approximation  was  suggested  (Rayleigh-aerosol  multiple  scattering 
ignored)  to  relate  water-leaving  radiance  (LK)  to  Lu  10 

LtW  =  UW  +  La  +  t(X)  LW(X)  +  LgUnt  (1  ) 

where  t(X)  is  diffuse  transmittance  of  the  atmosphere,  Ln  La  and  Lg\mt  are  radiance  contributions  due  to  Rayleigh,  aerosol, 
and  glint,  respectively.  Lr  is  derived  from  radiance  phase  functions  for  water  molecules  that  depend  on  incident  solar 
angles  (zenith  and  azimuth),  and  Fresnel  reflectance  estimates.  La  was  calculated  over  clcarwater  pixels  where  a 
minimum  water-leaving  radiance  at  665  nm  is  expected  (La  =  k  {Lt  (665)  -  Lr  (665)},  where  k  is  a  constant  assuming  a 
maritime  atmosphere) 10.  Lg/im  was  quantified  with  a  first-order  adjustment  by  subtracting  L{  (780)  to  Ltn.  Skylight  path 
radiance  contribution  was  assumed  small  due  to  the  relatively  thin  atmospheric  layer  between  the  airplane  and  the  sea 
surface.  Further  details  about  atmospheric  corrections  are  described  in  a  previous  work12.  Assuming  a  negligible 
attenuation  of  Ed  due  to  the  atmospheric  path  below  the  airplane,  the  remote  sensing  reflectance  above  the  sea-surface 
(Rrs(A,,  0f))  was  derived  as  normalized  water  leaving  radiance  (nZ,w(A,,  0+))  divided  by  E^XfiQ. 

In  case  1  waters,  R^,  0+)  can  be  approximately  related  (~20%  bias)  to  inherent  optical  properties  of  the  water  body 
with  the  following  expression1^: 


Rrs(^,  0+) «  0.54  R/Q  (2) 

R/Q  =  0.095  {bb  (X)/(bh  (X)  +  a  (X))}  (3) 

where  bb  is  the  total  backscattcring  coefficient  (water  +  particulates),  a  is  the  total  absorption  coefficient  of  water 
including  colored  dissolved  organic  matter  and  particulates  and  R/Q  is  a  shape  distribution  factor  that  is  influenced  by 
the  light  field  geometry. 

2.3  Active  optical  measurements 

Lidar  backscattering  measurements  were  obtained  with  a  Fish  Lidar  Oceanic  Experimental  (FLOE)  system14  mounted 
downwards  from  the  port  side  of  a  twin-engine  aircraft.  FLOE  was  set  up  15°off  vertical  to  minimize  specular  reflections 
from  the  sea.  For  each  lidar  pulse,  $  or  the  sum  of  water  (/(v)  and  particulate  (i.e.,  phytoplankton,  zooplankton,  fish)  (JJp) 
contributions  was  computed  from  photocathode  current  measurements  ( S)  as  a  function  of  depth  (z): 

S(z)  =  Afr(z)(L(zy2)e2az  +  B  (4) 

where  A  is  an  amplitude  parameter  that  depends  on  the  optical  system  parameters  and  the  geometry  (e.g.,  laser  pulse 
energy,  surface  losses,  receiver  area,  detector  responsivity),  L  is  the  optical  distance  from  the 

aircraft  to  the  measurement  depth  in  m,  a  is  the  lidar  attenuation  coefficient  in  m 4,  B  is  the  background  signal  level 
coming  from  skylight  contribution.  The  quantities  Aft*  and  a  were  found  for  each  lidar  pulse  based  on  equation  (4)  and 
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assuming  that  does  not  vary  with  depth,  and  ftp  is  zero  at  a  depth  of  2  m  and  at  the  maximum  penetration  depth  of 
each  lidar  pulse  (i.c.,  S(z)  above  10  standard  deviations  of  receiver  noise). 

For  most  of  the  surveyed  area,  FLOE  yielded  a  total  of  2,000  lidar  shots  (i.e.,  profiles)  per  minute  or  1  ‘scene’  in  about 
4. 1  km  distance.  This  corresponds  with  a  laser  sampling  rate  of  30  Hz  and  results  in  a  pixel  size  of  circa  2  m  (along- 
track)  by  5  m  (across-track)  by  0. 1  m  (along  the  vertical).  The  laser  is  linearly  polarized  and  has  beam  divergence  of  50 
mrad  during  daylight  hours  that  allow  P  measurements  as  deep  as  100  m.  However,  due  to  the  background  absorption  of 
laser  energy  with  depth,  FLOE  penetration  depth  in  our  study  area  was  30  m  depth  in  average  or  10*9  A  in  terms  of 
photocathode  current.  The  green  laser  source  was  pulsed  with  energy  of  100  mJ  and  a  length  of  15  nsec.  The  FLOE 
detector  has  a  cross-polarizer  in  front  of  the  telescope  to  maximizing  contrast  between  fish  and  smaller  light-scattering 
particles1  \ 

2.4  Calculation  of  remote  sensing  products 

To  examine  relationships  between  above-water  remote  sensing  reflectance  and  lidar  backscattering  measurements  we 
calculated  six  variables  based  on  MicroSAS  data  (5  spectral  band  ratios,  Rrs(410)/Rrs  (555),  Rrs(443)/Rrs(555), 
Rrs(490)/Rrs(555),  Rrs(508)/Rrs(555),  Rrs(443)/Rrs(490),  and  1  spectral  curvature  ratio,  G(l,l)  =  Rrs(490)2/{Rrs(443) 
Rre(508)}),  and  five  variables  based  on  FLOE-derived  /?  measurements  per  bin  (Maximum  P  per  bin,  /3Lax,  mean  /?  per 
bin,  /3tn,  standard  deviation  of  p  per  bin,  /4td,  integrated  p  per  bin,  $nt)  and  group  of  vertical  bins  (lidar  volume 
extinction  coefficient  of  P  between  0  and  5  m  depth,  ctp).  Notice  that  using  R^  ratios  allow  MicroSAS-FLOE 
comparisons  to  be  almost  independent  on  sunlight  illumination  conditions.  Likewise,  the  use  of  different  Rrs  wavelengths 
may  help  to  interpret  the  nature  of  optical  particulate  components  determining  /?(e. g.,  Rrs(443)/Rrs(490)  is  sensitive  to 
particle  size  distribution)16. 


2.5  Statistical  analysis 

Spatial  coherence  between  passive  and  active  optical  measurements  was  quantified  using  multiple  regression  analysis 
where  the  independent  variable  was  one  of  the  proposed  spectral  Rrs  ratios,  and  the  dependent  variables  were  the  /?- 
derived  parameters  calculated  at  different  depths. 

Before  each  run,  R^  and  p  parameters  were  aggregated  (MicroSAS  data  in  250-m  (along-track),  FLOE  data  in  250-m 
(along-track)  by  1  m  (along  the  vertical))  in  bins.  The  choice  of  250-m  spatial  resolution  in  the  horizontal  component 
coincides  with  the  maximum  spatial  resolution  provided  by  some  global  ocean  color  sensors  (c.g.,  MODIS).  Also,  this 
resolution  roughly  matches  the  swath  of  MicroSAS  radiometric  measurements  in  the  visible  range  (across-track  pixel 
size  =  200  m).  Intensity  of  covariation  between  MicroSAS  and  FLOE  variables  was  measured  based  on  semi-partial 
correlation  coefficients  (/>).  Relative  contribution  of  each  depth  to  spatial  changes  on  Rrb  ratios  was  estimated  by 
comparing  magnitude  of  p  of  only  those  depths  w  ith  significant  coefficients  at  95%  of  confidence  level. 


3.  RESULTS  AND  DISCUSSION 

The  simultaneous  use  of  passive  ocean  color  measurements  and  lidar  profiles  put  in  evidence  spatial  relationships 
between  Rrs  ratios  and  /^-derived  quantities  for  a  broad  range  of  spatial  scales  varying  between  250-m  and  50  km  along 
the  flight  direction  (Fig.  2,  Table  1). 
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Fig.  2  An  example  of  lidargram  during  August  17,  2002  survey  over  coastal  waters  of  Kodiak  Island.  The  whole  section 
encompasses  108  km  of  horizontal  distance  or  28  minutes  of  flight  duration.  Vertically  and  horizontally  integrated  lidar 
volume-backscattcring  (/?,„,)  is  plotted  in  log!0  scale  as  a  function  of  depth.  Spatial  integration  of  f3m  consisted  in  size 
bins  of  250  m  (horizontal)  by  1  meter  (vertical).  Rrs(490,  0f)/Rrs(553,  0+)  is  depicted  with  a  solid  line  and  relative  units. 
Remote  sensing  reflectance  ratio  spikes  in  along-track  bins  #  100,  1 10  and  250  were  related  to  patchy  distribution  of  low 
clouds  (i.e.,  drops  on  downwelling  irradiance  at  780  nm).  This  effect  was  amplified  in  those  R ratios  based  on  shorter 
wavelengths  in  the  numerator. 

In  Figure  2,  wc  show  the  best  spatial  coherence,  in  terms  of  p,  between  passive  and  active  optical  measurements 
performed  from  a  low  altitude  aircraft.  J3in(  and  Rrs(490,  0+)/Rrs(553,  0+)  values  were  positively  correlated,  and  that 
correlation  was  greater  in  the  first  meters  of  the  water  column  (Table  1).  Also,  spatial  variability  (magnitude)  of 
MicroSAS  Rrs  ratios  tended  to  be  higher  (lower)  for  locations  where  J3in(  decreased  drastically  with  depth  and  especially 
beyond  10  m  (e.g.,  bin  #  1  to  175).  Assuming  minor  surface  effects  (glint,  bubbles  and  foam),  the  intensification  of  lidar 
returns  near  the  sea  surface  was  likely  caused  by  greater  concentrations  of  phytoplankton  and  relatively  small  (<  250  m) 
fish  schools.  Surfacing  of  zooplankton  due  to  upward  migration  was  unlikely  during  the  sampling  period  since  these 
organisms  start  moving  toward  the  sea  surface  at  the  end  of  the  evening  twilight  (sunset  20:05  h,  local  time  in  Kodiak 
Island  during  August  17,  2002).  The  elevated  Rrs(490,  0+)/  Rrs(553,  0+)  associated  with  relatively  high  J3int  is  less 
straightforward  to  explain  due  to  the  lack  of  concurrent  in  situ  biological  data.  At  a  wavelength  of  532  nm,  zooplankton 
and  fish  have  a  stronger  interaction  with  the  lidar  waveform  in  terms  of  backscattcring  (i.e.,  more  reflective  targets) 
compared  to  phytoplankton  cells.  Likewise,  absorption  of  blue  wavelengths  (450-495  nm)  due  to  phytoplankton 
photosynthctic  pigments  and  colored  dissolved  organic  matter  is  a  major  light  attenuation  process  compared  those 
contributions  originated  from  fish  or  zooplankton.  Since  Rrs(553,  0+)  was  relatively  constant  between  bin  numbers  1-175 
and  175-345  (Fig.  3,  upper  curve),  the  main  change  on  Rrs(490,  0+)/  Rrs(553,  0+)  was  mainly  determined  by  blue  light 
attenuation  differences.  Given  that  fish  reflectivity  in  the  blue  is  relatively  low17,  the  increase  of  Rrs(490,  0+)  relative  to 
Rrs(553,  0+)  for  higher  J3int  values  was  probably  attributed  to  lower  phytoplankton  concentrations  (less  light  absorption  at 
490  nm)  and  greater  zooplankton  densities  (greater  retro-scattering  of  photons  at  532  nm).  Based  on  daily  FLOE 
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measurements,  integrated  ft  over  the  upper  20  m  and  averaged  over  a  distance  of  2.5  km  has  been  shown  to  be  positively 
related  to  zooplankton  settled  volume  during  the  same  month  of  200018. 

Table  1.  Statistical  relationships  between  above-water  remote  sensing  reflectance  ratios  and  under-water  lidar 
parameters  In  all  cases,  comparisons  between  optieal  measurements  were  performed  with  250-m  horizontal  resolution 
(MieroSAS  and  FLOE)  by  1  meter  vertical  resolution  (FLOE)  bins. 

flm\  arithmetic  average  of  lidar-derived  volume  baekseattering  (Jf)  per  bin  (sr'1  m'1),  pstd\  standard  deviation  of  /?per  bin 
(sr*1  m  ■'),  Pbu-  integrated  (3  per  bin  (m'1),  BinCorr:  lidar  depth  bins  showing  significant  semi-partial  correlation  with 
speetral  band  ratio  at  95%  eonfidenee  level,  BinMaxCorr:  lidar  depth  bin  having  maximum  semi-partial  correlation  with 
speetral  band  ratio,  ^ adf  adjusted  multiple  regression  coefficient. 


MieroSAS 
speetral  ratiod 

Lidar 

parameter 

BinCorr 

(m) 

BinMaxCorr 

(m) 

- 2 - 

r  adj 

1 

Psld 

1,2,7,8,12 

8 

0.19 

2 

Pm 

1,5,6,13 

1 

0.41 

Psld 

1,2,7,8,14 

1 

0.35 

Pin, 

1 

1 

0.44 

3 

Pm 

1 

1 

0.35 

Psld 

1,2,8,14 

1 

0.31 

Pin, 

1,3, 4, 5 

1 

0.41 

4 

Pm 

1 

1 

0.24 

Psld 

1,13 

1 

0.18 

Pin, 

1,20 

1 

0.26 

"Speetral  band  ratio  1 :  1^(443,  0+)/  Rrs(553,  0+),  2:  Rrs(490,  0+)/Rrs(553,  0+),  3:  Rrs(508,  0+)/Rrs(553,  0+),  and  4:  Rrs(443, 
0+)/  Rrs(490,  0(). 


Unlike  band  spectral  ratios,  spatial  variability  of  spectral  curvature  ratios  was  not  related  to  horizontal  variability  of  lidar 
parameters.  In  general,  G  (1,  1)  variations  reflect  change  on  scattering  versus  absorption  properties  of  phytoplankton 
communities  as  phytoplankton  blooms  develop12.  However,  the  weak  connection  between  G  and  lidar-baekscattcring 
parameters  was  more  likely  related  to  variability  of  at  443  nm  in  G  denominator  caused  by  ehanges  on  sea  surfaee 
roughness.  The  same  reason  may  explain  the  poorer  relationship  of  /3inh  j3m  and  /3s(d  with  Rrs(443,  0+)/  11^(553,  0f)  and 
Rrs(443,  0^)/Rrs(4()0,  0+)  (Table  1).  Lidar-based  parameters  sueh  as  (3max  were  also  highly  influenced  by  wind  and  wave 
effeets  near  the  sea  surfaee  becoming  a  variable  hard  to  prediet  based  on  above-water  water  leaving  radianee  ratios. 
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Fig.  3.  Spatial  variability  of  passive  and  active  remote  sensing  variables  that  link  above  and  in-water  optical  properties. 
In  left-y  axis  and  logarithm  scale  with  base  10,  Rrs(553,  0+):  above-water  remote  sensing  reflectance  derived  from 
MicroSAS  at  a  wavelength  of  553  nm  (sr  (upper  curve),  ap\  lidar  attenuation  coefficient  of  / 3 calculated  within  the  first 
five  meters  of  the  water  column  (m1)  (lower  curve). 

In  general  for  the  spectral  range  443-508  nm,  spatial  changes  on  Rrs  band  ratios  had  a  greater  association  with  pint  than 
with  pm ,  suggesting  a  primary  influence  of  number  of  particles  including  aggregates  with  respect  to  particle  composition 
modulating  horizontal  patterns  of  R^  and  /?  (Table  1,  Fig.  4A,C).  With  the  exception  of  Rrs(443,  0+)/Rrs(490,  0f) 
comparisons,  the  maximum  depth  at  which  lidar  parameters  still  have  significant  covariability  with  surface  Rrs  ratios  was 
commonly  greater  for  pstj  with  respect  to 

pmt  and  pm  (Fig.  4A,C-D).  This  should  not  be  surprising  since  pstd  had  generally  a  smaller  vertical  variation  compared  to 
Pmax,  Pint)  or  Pm >  thus  a  shorter  decorrelation  length  scale  as  a  function  of  depth  is  expected  between  MieroSAS-derived 
Rrs  ratios  and  lidar  parameters  related  with  abundance  or  type  of  optical  backscattcring  components  (Fig.  4D).  In 
general,  along-track  variability  of  Rrs  band  ratios  was  quite  indifferent  to  horizontal  changes  on  Cty(Fig.  3,  lower  curve). 
This  may  likely  related  to  a  weaker  influence  of  vertical  distribution  of  scatterers  compared  to  their  abundance  and  type 
on  Rrs  ratios  between  lidar  shots  and  for  p measured  within  0-5  depth. 
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Fig.  4.  Vertical  cross-section  of  lidar-derived  parameters  during  the  aerial  survey  made  in  August  17,  2002.  In  logl0 
scale,  A)  Vertically  integrated  volume  baekseattering,  fint  (sr1),  B)  Maximum  volume  backseattenng  per  bin,  (3max  (m*1 
sr*1),  C)  Mean  volume  baekseattering  per  bin,  (m  1  sr  !),  D)  Standard  deviation  of  volume  baekseattering  per  bin,  J3s((j 
(m  1  sr*1).  Each  bin  represents  250  m  along  the  flight  direction  and  1  m  section  as  a  function  of  water  depth. 


Interpretation  of  functionalities  between  Rrs  ratios  and  /?  spatial  variability  is  a  challenging  topic  due  mainly  to  two 
reasons.  Firstly,  even  using  the  same  wavelength,  the  origin  of  photons  measured  by  MieroSAS  (first  optical  depth)  and 
lidar  (0-30  m  depth  every  0. 1  m)  sensors  is  different,  and  secondly,  J3  is  more  influenced  by  relatively  ‘large’  targets  (e.g. 
zooplankton,  fish)  with  respect  to  Rrs  (e.g.,  phytoplankton)  since  lidar  is  a  more  collimated  light  source  with  respect  to 
the  sun.  Thus,  a  greater  fraction  of  target-reflectcd/medium-reflected  photons  occur  especially  when  particles  arc  larger. 
This  effect  is  expected  to  be  more  remarkable  at  shorter  wavelengths  (e.g  ,  <  450  nm)  and  may  explain  in  part  the 
absence  of  statistical  R^-/?  relationships  at  410  nm.  At  wavelengths  longer  than  the  FLOE  laser  wavelength,  the  Rrs-/? 
linkage  disappears  as  result  of  a  smaller  contribution  of  water-leaving  photons  in  the  red  spectrum  (i.e.,  >600  nm)  due  to 
the  light  absorption  by  water  near  the  sea  surface  (1  -2  m  depth). 


The  greatest  correlation  between  Rrs  ratios  and  (3  parameters  was  found  in  the  spectral  range  490-508  nm.  Within  this 
spectral  window,  fish  reflectivity  is  maximum  and  relatively  constant17,  thus  wc  suggest  that  observed  spectral 
differences  on  optical  relationships  were  mainly  accounted  by  particulate  components  other  than  fish  and  characterized 
by  having  a  major  influence  on  Rrs  (blue  light  absorption  by  phytoplankton)  and  (3  (green  light  baekseattering  by 
zooplankton)1*,  and  a  mutual  covariation  for  the  spatial  scale  under  study. 
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Our  preliminary  results  encourage  the  analysis  of  spatial  changes  on  Rrs  ratios  to  detect  sub-surface  optical  features  such 
that  drastic  change  on  /?  -parameters  at  12  m  depth  observed  after  bin  #  175.  The  use  of  Rrs  ratios  to  estimate  some 
aspects  of  vertical  distribution  of  lidar-derived  optical  properties  need  to  be  addressed  with  care  since  maximum  depth 
range  of  correlation  between  Rrs  ratios  and  fi  -parameters  varied  with  wavelength  and  lidar  backscattenng  property. 
Comparisons  between  shipboard  optical  profilers  (e.g.,  PRR),  airborne  Rrs,  and  lidar  back  scattering  parameters  at 
multiple  spatial  scales  (250  m  to  9  km)  arc  suggested  in  future  studies  to  better  understand  vertical  changes  in  optical 
properties  based  on  above-water  passive  optical  measurements. 
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